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Flat-A CDM Cosmogonies 
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Naoshi Sugiyama 7 , and Krzysztof M. Gorski 8 ' 9 

ABSTRACT 

We use Python I, II, and III cosmic microwave background anisotropy data 
to constrain cosmogonies. We account for the Python beamwidth and calibration 
uncertainties. We consider open and spatially- flat- A cold dark matter cosmogonies, 
with nonrelativistic-mass density parameter hi the range 0.1-1, baryonic-mass 
density parameter Qb in the range (0. 005-0. 029) h~ 2 , and age of the universe to in 
the range (10-20) Gyr. Marginalizing over all parameters but Qq, the combined 
Python data favors an open (spatially-flat-A) model with Qq ~ 0.2 (0.1). At the 2 a 
confidence level model normalizations deduced from the combined Python data are 
mostly consistent with those drawn from the DMR, UCSB South Pole 1994, ARGO, 
MAX 4 and 5, White Dish, and SuZIE data sets. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale 
structure of the universe 

1. Introduction 

Ganga et al. (1997a, hereafter GRGS) developed a technique to account for uncertainties, 
such as those in the beamwidth and the calibration, in likelihood analyses of cosmic microwave 

department of Physics, Kansas State University, Manhattan, KS 66506. 

2 Centro de Astrofi'sica da Universidade do Porto, Rua das Estrelas s/n, 4100 Porto, Portugal. 
3 Center for Particle Astrophysics, University of California, Berkeley, CA 94720. 
4 Copernicus Astronomical Center, Bartycka 18, 00-716 Warszawa, Poland. 

Infrared Processing and Analysis Center, MS 100-22, California Institute of Technology, Pasadena, CA 91125. 
6 Snow and Ice Research Group, University of Nebraska, Lincoln, NE 68583-0850. 

7 Department of Physics, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-ku, Kyoto 606-8502, Japan. 
8 Theoretical Astrophysics Center, Juliane Maries Vej 30, 2100 Copenhagen 0, Denmark. 
9 Warsaw University Observatory, Aleje Ujazdowskie 4, 00-478 Warszawa, Poland. 



- 2 - 



background (CMB) anisotropy data. This technique has been used in conjunction with 
theoretically-predicted CMB anisotropy spectra in analyses of the Gundersen et al. (1995) UCSB 
South Pole 1994 data, the Church et al. (1997) SuZIE data, the MAX 4+5 data (Tanaka et al. 
1996; Lim et al. 1996), the Tucker et al. (1993) White Dish data, and the de Bernardis et al. 
(1994) ARGO data (GRGS; Ganga et al. 1997b, 1998; Ratra et al. 1998, 1999a, hereafter R99a). 
A combined analysis of all these data sets is presented in Ratra et al. (1999b, hereafter R99b). 

In this paper we present a similar analysis of CMB anisotropy data from the Python I, II, 
and III observations performed at the South Pole (Dragovan et al. 1994, hereafter D94; Ruhl et 
al. 1995b, hereafter R95; Piatt et al. 1997, hereafter P97). The Python detectors and telescope 
are described by Ruhl (1993) and D94; also see Ruhl et al. (1995a) and Alvarez (1996). In what 
follows we review the information needed for our analysis. 

Python I, II, and III CMB data were taken in a frequency band centered at 90 GHz with four 
bolometric detectors centered at the corners of a 2? 75 by 2? 75 square on the sky. The beam profiles 
are well-approximated by a Gaussian of FWHM 0?75 ± 0?05 (one standard deviation uncertainty). 
Observations were centered at a = 23 1 . 1 37, 5 = —49? 44 (J2000.0). Python I and II data were taken 
at a single telescope elevation. Python III data were taken at this fiducial elevation as well as two 
additional elevations offset 2?75/3 on the sky above and below the fiducial elevation. The reduced 
Python data are shown in Figure 1. 

All of the Python measurements were made by switching the four beams horizontally across 
the sky in a three-point pattern by rotating a vertical flat mirror at 2.5 Hz. This chopping pattern 
was then combined with slow (typically 0.1 Hz) azimuthal beam switching of the entire telescope 
to produce a four-beam response to a sky signal. 

The chopper throw and azimuthal telescope beam switching were both 2?75 on the sky for 
the Python I and II observations. Python I (hereafter I) resulted in 16 data points, 8 each from 
the lower and upper rows of detectors (D94). Python II observed two sets of points on the sky 
(R95), although one of the detectors did not work during this season. The first set of observations, 
hereafter IIA, overlapped the I points and yielded 7 measurements from the lower row of detectors 
and 8 from the upper row. The second set, hereafter IIB, measured the same number of points on 
the sky as the first set, but these were offset in azimuth relative to the I points by — 2?75/2 on the 
sky. 

Two series of measurements were made at each of the three elevations observed during the 
Python III season (P97). For each series, the physical throw of the chopper was the same at all 
elevations. Because the actual throw on the sky depends on elevation, the Python III beam was 
smeared to a Gaussian FWHM of 0?82 ± 0?05 (one standard deviation uncertainty) when forming 
the four-beam pattern to account for the imperfect overlap of the beams caused by this effect. 
This procedure also accounted for the relative pointing uncertainty and fluctuations in the chopper 
throw. 

The first series of Python III measurements, hereafter IIIL, used the same chopper and beam 
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switch parameters as Python I and II. The IIIL data consists of two sets of points taken at the 
I elevation but offset in right ascension by —2? 75/3 and —2 x 2? 75/3 on the sky. Two sets of 
measurements were also made at each of the lower and upper elevations, but these were offset in 
right ascension by 0?0 and +2 x 2?75/3 on the sky relative to I. 

The second series of Python III measurements, hereafter HIS, were made with both the 
chopper throw and telescope beam switch reduced to 2?75/3. For each of the three telescope 
elevations, the HIS data consists of points separated horizontally by 2?75/3 on the sky for each of 
the 2 rows of detectors. Figure 1 shows the points observed by IIIL and HIS at each elevation. 
Together, I, II, and III densely sample a 5?5 by 22° region of the sky. 

The 1 a absolute calibration uncertainty in the Python data is 20% (D94; R95; P97) and is 
accounted for in our analysis. The absolute pointing uncertainty is 0?1 (D94; R95) and is not 
accounted for in our analysis. 

In §2 we summarize the computational techniques used in our analysis. See GRGS and R99a 
for detailed discussions. Results are presented and discussed in §3. Conclusions are given in §4. 

2. Summary of Computation 

The zero-lag window function for the Python observations are shown in Figure 2 and the 
zero-lag window function parameters are in Table 1. 

In this paper we focus on a spatially open CDM model and a spatially flat CDM model 
with a cosmological constant A. These low density models are largely consistent with current 
observational constraints.^ For recent discussions see Park et al. (1998), Retzlaff et al. (1998), 
Croft et al. (1999), and Peebles (1999). 

The models have Gaussian, adiabatic primordial energy-density power spectra. The 
fiat-A model CMB anisotropy computations use a scale-invariant energy-density perturbation 
power spectrum (Harrison 1970; Peebles & Yu 1970; Zel'dovich 1972), as predicted in the 
simplest spatially-flat inflation models (Guth 1981; Kazanas 1980; Sato 1981). The open model 
computations use the energy-density power spectrum (Ratra & Peebles 1994, 1995; Bucher, 
Goldhaber, & Turok 1995; Yamamoto, Sasaki, &; Tanaka 1995) predicted in the simplest 
open-bubble inflation models (Gott 1982). The computation of the CMB anisotropy spectra is 
described by Stompor (1994) and Sugiyama (1995). 

As discussed in R99a, the spectra are parameterized by their quadrupole-moment amplitude 



1 While not considered in this paper, a time-variable cosmological "constant" dominated spatially-flat model is 
also largely consistent with current data (e.g., Peebles & Ratra 1988; Sugiyama & Sato 1992; Ratra & Quillen 1992; 
Frieman & Waga 1998; Ferreira & Joyce 1998; Wang & Steinhardt 1998; Carroll 1998; Hu et al. 1999; Huterer & 
Turner 1999; Liddle & Scherrer 1999; Starobinsky 1998). 
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Qrms-PSj the nonrelativistic-mass density parameter Qq, the baryonic-mass density parameter Qb, 
and the age of the universe to- The spectra are computed for a range of 0,q spanning the interval 
0.1 to 1 in steps of 0.1, for a range of fish 2 [the Hubble parameter h = ifo/(100 km s _1 Mpc -1 )] 
spanning the interval 0.005 to 0.029 in steps of 0.004, and for a range of to spanning the 
interval 10 to 20 Gyr in steps of 2 Gyr. In total 798 spectra were computed to cover the 
cosmological-parameter spaces of the open and flat- A models. Figure 2 shows examples of the 
CMB anisotropy spectra used in our analysis. Other examples are in Figure 2 of R99a and Figure 
1 of R99b. 

While it is of interest to also consider other cosmological parameters, such as tilt or gravity 
wave fraction or a time- variable cosmological "constant" (instead of a constant A), to make the 
problem tractable we have focussed on the four parameters mentioned above. We emphasize 
however that the results of the analysis are model dependent. For instance, a time-variable A 
model would likely lead to a different constraint on Qq than that derived below in the constant A 
model. 

Following GRGS, for each of the 798 spectra considered the "bare" likelihood function is 
computed at the nominal beamwidth and calibration, as well as at a number of other values of 
the beamwidth and calibration determined from the measurement uncertainties. The likelihood 
function used in the derivation of the central values and limits is determined by integrating 
(marginalizing) the bare likelihood function over the beamwidth and calibration uncertainties 
with weights determined by the measured probability distribution functions of the beamwidth 
and the calibration. See GRGS for a more detailed discussion. When marginalizing over the 
beamwidth uncertainty we have checked in a few selected cases that the five-point Gauss-Hermite 
quadrature summation approximation to the integral agrees extremely well with the three-point 
Gauss-Hermite approximation used by GRGS (and for most of the analysis in this paper). The 
likelihoods are a function of four parameters mentioned above: Qrms— PS, ri , ^Bh 2 , and to- We 
also compute marginalized likelihood functions by integrating over one or more of these parameters 
after assuming a uniform prior in the relevant parameters. The prior is set to zero outside the 
ranges considered for the parameters. GRGS and R99a describe the prescription used to determine 
central values and limits from the likelihood functions. In what follows we consider 1, 2, and 3 a 
highest posterior density limits which include 68.3, 95.4, and 99.7% of the area. 



3. Results and Discussion 

Table 2 lists the derived values of Qrms-PS an d bandtemperature bT\ for the flat bandpower 
spectrum, for various combinations of the I, II, and III data. These numerical values account 
for the beamwidth and calibration uncertainties. The last two STi entries in Table 2 are quite 
consistent with those derived by P97; the small differences reflect the different methods used to 
account for beamwidth and calibration uncertainties here and in P97. 



- 5 - 



For the flat bandpower spectrum the combined I, II, and III data average 1 a 5T[ error bar 
is ~ 25%0 : Python data results in a very significant detection of CMB anisotropy, even after 
accounting for beamwidth and calibration uncertainties. Note that the calibration uncertainty, 
20%, is the most important contributor to this error bar. 

A number of other interesting conclusions follow from the entries in Table 2. Comparing the 
I+II and IIIL results, which are from experiments which probe almost identical angular scales, we 
see that the IIIL amplitude is ~ 1 a higher than the I+II amplitude. Comparing the result from 
the analysis of the coadded I and IIA data (which are from experiments with identical window 
functions) and the result from the full analysis of the I and IIA data (which takes into account all 
the spatial correlations), we see that the deduced amplitudes are almost identical. This is probably 
mostly a reflection of the fact that the individual I and IIA amplitudes are almost identical (see 
Table 2). 

As discussed in R99a and R99b, the four-dimensional posterior probability density distribution 
function L((3rms-PS> ^0; &Bh 2 , to) is nicely peaked in the Q r ms-PS direction but fairly flat in 
the other three directions. Marginalizing over Qi-ms-PS results in a three-dimensional posterior 
distribution L(Qq, Qgh 2 , to) which is steeper, but still relatively flat. As a consequence, limits 
derived from the four- and three-dimensional posterior distributions are generally not highly 
statistically significant. We therefore do not show contour plots of these functions here. 
Marginalizing over Qrms-PS and one other parameter results in two-dimensional posterior 
probability distributions which are more peaked. See Figures 3 and 4. As in the ARGO (R99a) 
and combination (R99b) data set analyses, in some cases these peaks are at an edge of the 
parameter range considered. 

Figure 3 shows that the two-dimensional posterior distributions allow one to distinguish 
between different regions of parameter space at a fairly high formal level of confidence. For 
instance, the open model near S7o ~ 0.75, f2#/i 2 ~ 0.03, and to ~ 20 Gyr, and the flat-A model 
near f^o ~ 0.6, f2#/i 2 ~ 0.03, and to ~ 20 Gyr, are both formally ruled out at ~ 3 a confidence. 
However, we emphasize, as discussed in R99a and R99b, care must be exercised when interpreting 
the discriminative power of these formal limits, since they depend sensitively on the fact that the 
uniform prior has been set to zero outside the range of the parameter space we have considered. 

Figure 4 shows the contours of the two-dimensional posterior distribution for Q r ms-PS an d 
Qo, derived by marginalizing the four-dimensional distribution over f2#/i 2 and to- These are shown 
for the combined Python I, II, and III data, the DMR data, and three combinations of data from 
the SP94, ARGO, MAX 4+5, White Dish, and SuZIE experiments (R99b), for both the open and 
flat-A models. Constraints on these parameters from the combined Python data are consistent 
with those from the DMR data for the flat-A models, panel a), while for the open model, panel b), 



11 For comparison, the corresponding DMR 1 a STi error bar is ~ 10 — 12% (depending on model, Gorski et al. 
1998). 
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consistency at 2 a (1 a) requires fio ~ 0.2 (0.35). The combined Python data amplitudes are a 
little higher than those derived from the other small-scale data combinations, panels c) — h), but 
at 2 a confidence the various amplitudes are mostly consistent. 

Figure 5 shows the one-dimensional posterior distribution functions for Qq, £l B h 2 , to, and 
Qrms-PS> derived by marginalizing the four-dimensional posterior distribution over the other three 
parameters. From these one-dimensional distributions, the combined I, II, and III data favors an 
open (flat-A) model with Q = 0.19 (0.10), or n B h 2 = 0.005 (0.005), or t = 10 (13) Gyr, amongst 
the models considered. At 2 a confidence the combined Python data formally rule out only small 
regions of parameter space. From the one-dimensional distributions of Figure 5, the data requires 
n < 0.71 or > 0.8 (Qo < 0.55 or > 0.63), or n B h 2 < 0.028 (n B h 2 < 0.028), or to < 20 Gyr (t 
< 20 Gyr) for the open (flat-A) model at 2 a. As discussed in R99a and R99b, care is needed 
when interpreting the discriminative power of these formal limits. These papers also discuss a 
more conservative Gaussian posterior distribution limit prescription. Using this more conservative 
prescription, we find only an upper 1 a limit on Oo ( ~ 0.5) in the open model. 

While the statistical significance of the constraints on cosmological parameters is not high, 
it is reassuring that the combined Python data favor low-density, low Q B h 2 , young models, 
consistent with some of the indications from the combinations of CMB anisotropy data considered 
by R99b, and the indications from most recent non-CMB observations (see discussion in R99b). 

The peak values of the one-dimensional posterior distributions shown in Figure 5 are listed 
in the figure caption for the case when the four-dimensional posterior distributions are normalized 
such that L(Q rms _ps = fiK) = 1. With this normalization, marginalizing over the remaining 
parameter the fully marginalized posterior distributions are 2 x 10 106 (1 x 10 106 ) for the open 
(flat-A) model and the combined Python data. This is qualitatively consistent with the indication 
from panels a) and b) of Figure 5 that the most-favored open model is somewhat more favored 
than the most-favored flat-A one. 



4. Conclusion 

The combined Python I, II, and III data results derived here are mostly consistent with 
those derived from the DMR, SP94, ARGO, MAX 4+5, White Dish and SuZIE data. The 
combined Python data significantly constrains Qrms-PS (f° r the flat bandpower spectrum 
Qrms-PS = 40 jiK at 1 a) and weakly favors low-density, low Q B h 2 , young models. 

We acknowledge helpful discussions with D. Alvarez, M. Dragovan, G. Griffin, J. Kovac, and 
J. Ruhl. This work was partially carried out at the Infrared Processing and Analysis Center 
and the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the 
National Aeronautics and Space Administration. KG also acknowledges support from NASA ADP 
grant NASA-1260. BR and GR acknowledge support from NSF grant EPS-9550487 with matching 
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AISRP grant NAG-3941 and help from Polish Scientific Committee (KBN) grant 2P03D00813. 
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Table 1: Numerical Values for the Zero-Lag Window Function Parameters a 





l e -0.5 


l c 


l m i e -0.6 




Python I/II 


53 


91.7 


73 99 


1.34 


Python IIIL 


52 


87.7 


72 98 


1.30 


Python HIS 


128 


171 


176 230 


0.623 



a The value of I where Wi is largest, i m , the two values of I where Wi _ a s = e 0,5 Wi m , i e -o.5, the effective multipole, 
Z c = /(IWi)//(Wi). and I(Wi) = Ei= 2 (' + 0.5)W,/{J(J + 1)}. 
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Table 2: Numerical Values for Qrms-PS an d from Likelihood Analyses Assuming a Flat 
Bandpower Spectrum 



Data Set a 


Qrms— PS 

(11K) 


Ave. Abs. Err. c 
(uK) 


Ave. Frac. Err. d 


m h 

(MK) 




LR C 


I 


ay 28 


14 


36% 


61 


87 
43 


2 


x 10 14 


IIA 


39 27 


15 


38% 


60 


no 

43 


9 


x 10 12 


C(I+IIA) 


40 57 

u 29 


14 


35% 


62 


89 
45 


3 


x 10 31 


I+IIA 


40 56 

w 29 


14 


35% 


61 


87 
45 


2 


x 10 31 


IIB f 


9.3 2 ° 


10 


110% 


14 


31 





1 


II 


01 23 


10 


34% 


48 


68 
35 


1 


x 10 12 


C(I+IIA)+IIB 


34 47 

°^ 26 


11 


31% 


53 


73 
40 


1 


x 10 30 


I+II 


34 46 

°^ 25 


11 


31% 


52 


72 
39 


1 


x 10 30 


IIIL 


41 55 

* 1 32 


11 


28% 


64 


84 
49 


3 


x 10 33 


HIS 


42 56 

4Z 33 


11 


27% 


66 


87 
51 


2 


x 10 29 


III 


41 54 

^ x 33 


10 


25% 


64 


83 
51 


6 


x 10 68 


I+II+IIIL 


3Q 51 
31 


10 


27% 


60 


80 
47 


3 


x 10 67 


I+II+III 


40 52 
w 32 


10 


25% 


63 


81 
50 


1 x 10 105 



a C(.. .+...) means that the data from the two sets of observations have been coadded prior to analysis, II refers to the 

combined IIA and IIB data, and III refers to the combined IIIL and HIS data. 
b For each data set, the first of the three entries is where the posterior probability density distribution function peaks 

and the vertical pair of numbers are the ±1 a (68.3% highest posterior density) values. 
c Average absolute error on Q r ms-ps in fjK- 
d Average fractional error, as a fraction of the central value. 
e Likelihood ratio. 

f IIB does not have a 2 cr highest posterior density detection; the appropriate equal tail 2 a upper limits are 50 [iK 
(Q r m S -ps) and 77 fiK (6T t ). 
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Figure Captions 

Fig. 1.- Measured thermodynamic temperature differences (with ±l-cx error bars). Panel a) shows 
Python I (open circles), IIA (filled squares, offset horizontally from true positions for clarity), IIB 
(open squares), and IIIL (crosses) data, while panel b) shows Python HIS data (crosses). 

Fig. 2.- CMB anisotropy multipole moments 1(1 + l)Q/(27r) x 10 10 (solid lines, scale on left 
axis, note that these are fractional anisotropy moments and thus dimensionless) as a function of 
multipole /, for selected models normalized to the DMR maps (Gorski et al. 1998; Stompor 1997). 
Panels a) — c) show selected flat- A models. The heavy lines are the VLq = 0.1, Qsh 2 = 0.005, and 
to = 12 Gyr case, which is close to where the combined Python data likelihoods (marginalized over 
all but one parameter at a time) are at a maximum. Panel a) shows five fig/i 2 = 0.005, to = 12 
Gyr models with Qq = 0.1, 0.3, 0.5, 0.7, and 0.9 in descending order at the I ~ 200 peaks. Panel 
b) shows seven fl = 0.1, t = 12 Gyr models with Q, B h 2 = 0.029, 0.025, 0.021, 0.017, 0.013, 0.009, 
and 0.005 in descending order at the I ~ 200 peaks. Panel c) shows six Qq = 0.1, VLsh? = 0.005 
models with to = 20, 18, 16, 14, 12, and 10 Gyr in descending order at the I ~ 200 peaks. Panels 
d) — f) show selected open models. The heavy lines are the ^o = 0.2, Vtsh 2 = 0.005, and to = 10 
Gyr case, which is close to where the combined Python data likelihoods (marginalized over all but 
one parameter at a time) are at a maximum. Panel d) shows five fi^/i 2 = 0.005, to = 10 Gyr 
models with fio = 1> 0.8, 0.6, 0.4, and 0.2 from left to right at the peaks (the peak of the = 0.2 
model is off scale). Panel e) shows seven $7o = 0.2, to = 10 Gyr models with Og/i 2 = 0.029, 0.025, 
0.021, 0.017, 0.013, 0.009, and 0.005 in descending order at I ~ 400. Panel /) shows six ft = 0.2, 
VL B h 2 = 0.005 models with t = 20, 18, 16, 14, 12, and 10 Gyr in descending order at I ~ 400. 
Also shown are the Python zero-lag window functions W\ (scale on right axis): I/II (long-dashed 
lines), IIIL (short-dashed lines), and HIS (dotted lines). See Table 1 for W;-parameter values. 

Fig. 3.- Confidence contours and maxima of the combined Python data two-dimensional posterior 
probability density distribution functions, as a function of the two parameters on the axes of each 
panel (derived by marginalizing the four-dimensional posterior distributions over the other two 
parameters). Dashed lines (crosses) show the contours (maxima) of the open case and solid lines 
(solid circles) show those of the flat-A model. Panel a) shows the (f2#/i 2 , Qq) plane, and panel b) 
shows the (to, fio) plane. 

Fig. 4.- Confidence contours and maxima of the two-dimensional (Qrms-PS) ^o) posterior 
probability density distribution functions. Panels a), c), e), & g) in the left column show the flat-A 
model and panels 6), d), f), & h) in the right column show the open model. Note the different scale 
on the vertical ((Jrms-ps) axes of pairs of panels in each row. Heavy lines show the ±1 and ±2 a 
confidence limits and solid circles show the maxima of the two-dimensional posterior distributions 
derived from the combined Python I, II, and III data. Shaded regions show the two-dimensional 
posterior distribution 1 a (denser shading) and 2 a (less dense shading) confidence regions for the 
DMR data (Gorski et al. 1998; Stompor 1997) in panels a) & b); for the SP94, ARGO, MAX 4 
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and 5, White Dish and SuZIE data combination (R99b) in panels c) Sz d); for the previous data 
combination excluding SuZIE (R99b) in panels e) & /); and for the SP94Ka, MAX 4 ID, and 
MAX 5 HR data combination (R99b) in panels g) & h). The DMR results are a composite of 
those from analyses of the two extreme data sets: i) galactic frame with quadrupole included and 
correcting for faint high-latitude galactic emission; and ii) ecliptic frame with quadrupole excluded 
and no other galactic emission correction (Gorski et al. 1998). In panels c) — h) crosses show the 
maxima of the appropriate non-Python data two-dimensional posterior distributions. 

Fig. 5.- One-dimensional posterior probability density distribution functions for f^o, ^Bh 2 , to, 
and Q r ms-PS (derived by marginalizing the four-dimensional one over the other three parameters) 
in the open and flat-A models. These have been renormalized to unity at the peaks. Dotted 
vertical lines show the confidence limits derived from these one-dimensional posterior distributions 
and solid vertical lines in panels g) and h) show the ±1 and ±2 a confidence limits derived by 
projecting the combined Python I, II, and III data four-dimensional posterior distributions. The 
2 a DMR (marginalized and projected) confidence limits in panels g) and h) are a composite of 
those from the two extreme DMR data sets (see caption of Figure 4). When the four-dimensional 
posterior distributions are normalized such that L(Q rm s-PS = fiK) = 1, the peak values of the 
one-dimensional distributions shown in panels a) — h) are 2 x 10 106 , 3 x 10 106 , 6 x 10 107 , 7 x 10 107 , 
1 x 10 105 , 2 x 10 105 , 9 x 10 104 , and 6 x 10 104 , respectively. 
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